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ABSTRACT 

We analyze the correlation of QCD one-loop effects on the partial widths of 
the three neutral Higgs bosons of the MSSM decaying into quark-antiquark 
pairs. The SUSY-QCD effects turn out to be comparable or even larger than 
the standard QCD effects and are slowly decoupling in a wide window of the 
parameter space. Our study is aimed at elucidating the possible supersymme- 
tric nature of the neutral Higgs bosons that might be discovered in the near 
future at the Tevatron and/or at the LHC. In particular, we point out the 
presence of potentially large SUSY corrections to the various neutral Higgs 
production cross-sections. 



To glimpse into the relevance of addressing the issue of the width of a Higgs boson, 
notice that if a heavy neutral Higgs is discovered and is found to have a narrow width, it 
would certainly not be the SM Higgs, whilst it could be a SUSY Higgs. In fact, a heavy 
enough SM Higgs boson is expected to rapidly develop a broad width through decays into 
gauge boson pairs whereas the SUSY Higgs bosons cannot in general be that broad since 
their couplings to gauge bosons are well-known to be suppressed [jlj . In compensation, 
their couplings to fermions (especially to heavy quarks) can be considerably augmented. 
Thus the width of a SUSY Higgs should to a great extent be given by its hadronic width; 
even so a heavy H° and A is in general narrower than a SM Higgs of the same mass. 
Alternatively, if the discovered neutral Higgs is sufficiently light that it cannot decay into 
gauge boson pairs, its decay width into relatively heavy fermion pairs such as r + r _ , and 
especially into bb, could be much larger than that of the SM Higgs, because of tan/3- 
enhancement of the fermion couplings |lj . Hence, it becomes clear that the hadronic 
width may play a very important role in the study of the MSSM higgses, already at the 
tree-level. 

As is well-known [p]] , a generic two-doublet Higgs sector is composed of a charged Higgs 
pseudoscalar boson, H ± 1 a neutral CP-odd ( "pseudoscalar" ) boson, A , and two neutral 
CP-even Higgs "scalars", h° and H° (M h o < M H o). In the specific case of the MSSM, 
and because of the supersymmetric constraints, only two parameters, e.g. (M^o, tan/3), 
are independent at the tree-level Therefore, in the MSSM definite predictions can be 
made and tested. 

The aim of this work is to complete the analysis of the strong SUSY corrections to 
the hadronic decay widths of the Higgs bosons of the MSSM that we have initiated in the 
companion paper 0. In the latter reference whose notation and definitions we assume 
hereafter we have defined two interesting domains of the general MSSM parameter, the 
so-called Regions I and II, where the physics of the supersymmetric Higgs bosons can be 
especially relevant. Depending on the region of parameter space considered, not all the 
Higgs particles of the MSSM are allowed to decay hadronically in a significant manner. 
On the one hand, the process H + — > tb, which requires Mh+ > m t + mj, ~ 180 GeV, is 
permitted in Region II and the SUSY-QCD corrections can be relevant in that region 0. 
When H + — > tb is allowed, such a decay is by far the main hadronic decay mode of a 
SUSY charged Higgs boson. If, however, the condition M H + > m t + mb is not satisfied, the 
remaining hadronic decays available to H + are not so appealing since the corresponding 
branching fractions lie below the leptonic r-mode H + — > r + z/ T . On the other hand, we 
may turn our attention to the various hadronic neutral Higgs decays $ l — > qq ($ 1 = 
A , $ 2 = h°, $ 3 = H°). Of these, we will neglect the decays leading to light qq final 
states since their branching ratio is very small. Thus, for the lightest neutral Higgs, h°, we 
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will concentrate on just the decay h° — > bb, whereas for A , H° (which can be arbitrarily 
heavy) we shall consider the channels A , H° — > bb and A , H° —* ti. 

Moreover, should the physical domain of the MSSM parameter space turn out to lie 
in Regions I or II, then we shall see that the hadronic widths of the MSSM Higgs bosons 
must incorporate important virtual SUSY signatures. The latters could be extracted from 
measured quantities by subtracting the corresponding conventional QCD effects, which 
can be easily computed by adapting the results of Refs.|3|, |J. Knowledge of the SUSY 
corrections could be determinant to trace the nature of those scalars and establish whether 
they are truly supersymmetric Higgs particles. 

Although we elaborate here mostly on the Higgs strategies at hadron colliders, such 
as the Tevatron and especially the LHC, it should be clear that the kind of effects that we 
wish to study have an impact on Higgs physics in e + e~ machines as well, where the neutral 
Higgs states can be produced through e.g. e + e~ — > Z h°(H°) and e + e~ — > A h°(H°). The 
observed cross-sections for these processes are equal to the production cross-sections times 
the Higgs branching ratios. Thus, in an e + e~ environment one aims more at a measure- 
ment of the various branching ratios (or, more precisely: ratios of branching ratios) of the 
fermionic Higgs decay modes rather than of the partial widths themselves. For instance, in 
e + e~ one would naturally address the measurement of BR(<& 1 — > bb) / BR(<& % — ► t + t~); 
in fact, this observable should receive large SUSY-QCD corrections if — > bb proves to 
be, as we shall see, very sensitive to the strong supersymmetric effects. 

In hadron machines an actual measurement of the hadronic partial widths and in 
general of the effective hadronic vertices <3> l qq(q = t, b) should be feasible. Let us briefly 
remind of the five basic mechanisms for neutral Higgs production in a hadron collider |5j . 
They have been primarily described for the SM Higgs, Hg M , but can be straightforwardly 
extended to the three neutral higgses, $ l , of any two-doublet Higgs sector (see Fig. la. for 
a sketch of some of these mechanisms): 

• (i) Gluon-gluon fusion: g g — > $ l ; 

• (ii) WW(ZZ) fusion: qq-> qq&; 

• (iii) Associated W(Z) production: qq — > W(Z) 

• (iv) ti fusion: g g — > and 

• (v) bb fusion: gg — > 66$*. 

It has been known for a long time that in the SM, where only one neutral Higgs Hg M 
is present, mechanism (i) provides the dominant contribution over most of the accessible 
range. For very large (obese) SM Higgs mass, however, mechanism (ii) eventually takes 
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over; the rest of the mechanisms are subleading, and in particular bb fusion is negligible 
in the SM. 

Remarkably enough, this situation could drastically change in the MSSM. As noted 
above for the Higgs decays, also the production mechanisms of the MSSM Higgs scalars 
can be very different from the SM [H. For instance, whereas one-loop g g-fusion in the 
SM is dominated by a top quark in the loop, this is not always so in the MSSM where the 
new couplings turn out to enhance, at high tan /3, the 6-quark loops and make them fully 
competitive with the top quark loops (Fig. la). Furthermore, mechanism (ii) becomes 
suppressed by the SUSY couplings; e.g. in Region I the lightest neutral Higgs, h°, has 
a very small coupling to the weak gauge bosons as compared to H SM . In this respect 
the situation with the CP-odd scalar, A , is even worse, for it can never be produced by 
mechanism (ii) at the tree-level. In contrast, bb fusion (Fig. la), which was negligible in 
the SM, can be very important in the MSSM at large tan (3, especially in Region I but 
also in Region II. As a matter of fact, for large enough tan/?, the 66-fusion cross-section 
can be larger than that for any mechanism for producing a SM Higgs boson of similar 
mass . 

Our interest in the production mechanisms mentioned above stems from the fact that 
the radiative effects could play a crucial role. This is true already in the SM. For example, 
the conventional QCD corrections to g g — > H SM , which is the dominant process for the 
production of a light and an intermediately heavy Higgs boson, are known to be large 0. 
A similar conclusion holds for an obese SM Higgs boson produced at very high energies by 
means of the WW(ZZ) -fusion mechanisms; here, again, non-negligible radiative effects 
do appear ||. Therefore, the production cross-section for H SM is expected to acquire 
valuable quantum corrections both for light and for heavy Higgs masses. This is not so, 
however, for the corresponding width. In fact, only for a heavy SM Higgs, namely, with 
a mass above the vector boson thresholds, the corrections to its decay width can be of 
interest; for a light SM Higgs, instead, light enough that it cannot decay into gauge boson 
pairs, the decay width is very small and thus the corresponding quantum effects are of no 
practical interest. 

Now, in contradistinction to the SM case, the hadronic vertices <3> J qq(q = t, b) could be 
the most significant interactions for MSSM higgses, irrespective of the value of the Higgs 
masses. In fact, these vertices can be greatly enhanced. Therefore, if large radiative 
corrections may modify the effective structure of these interactions, it is clear that they 
should be taken into account and could be of much practical interest. In what follows we 
shall substantiate that in the MSSM the bb and $ l tt vertices involved in mechanisms 
(i), (iv) and (v) above could receive very large SUSY-QCD corrections (in some cases 
above 50%) and so, if these effects are there, they will be reflected in the Higgs boson 
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partial widths and to a large extent also in the production cross-sections. In this respect 
the aforementioned b 6-fusion mechanism, which is highly operative at large tan /3, can be 
very sensitive to these SUSY-QCD corrections. To our knowledge, these matters have 
not been discussed in the literature and could play a momentous role to decide whether 
a neutral Higgs hypothetically produced in a hadron collider is super symmetric or not. 

While it goes beyond the scope of this note to compute the SUSY-QCD corrections to 
the production processes themselves, we have performed a detailed analysis of the partial 
decay widths, which are the canonical observables that should be first addressed to probe 
the new quantum corrections to the basic interaction vertices. In this way, a definite 
prediction is made on the properties of a physical observable, and moreover this should 
suffice both to exhibit the relevance of the SUSY quantum effects and to demonstrate the 
necessity to incorporate these corrections in a future, truly comprehensive, analysis of the 
cross-sections, namely, an analysis where one would include the quantum effects on all 
the relevant production mechanisms within the framework of the MSSM. 

Let us now concentrate on the diagrams depicted in Fig. lb. Since we adopt the 
same framework as in the companion paper ||, we shall obviate all the unessential details 
already defined there. The interaction Lagrangian describing the $ J q q- vertex in the 
MSSM is: 

U qq = [<(Q)PL + a l R {q)P R ] q. (1) 

We shall focus on top and bottom quarks (q = t,b). In a condensed and self-explaining 
notation we have defined 

a R (t,b) = — a^(£, b) = (i cot /3, i tan/3) , 
a 2 R (t,b) = a 2 L (t,b) = (-c a /sj3,s a /cp) , 

a 3 R (t,b) = a 3 L (t,b) = (s a /s p ,-c a /cp) , (2) 

with c a = cos a, s@ = sin (3 etc. (Angles a and j3 are related in the usual manner prescribed 
by the MSSM |T[].) Apart from the SUSY-QCD interactions involving gluinos and squarks, 
a very relevant piece of our calculation is the interaction Lagrangian between neutral 
higgses and squarks. In compact notation, it can be cast as follows: 

^=^* i %<#l>*, (3) 
where we have introduced the mass-eigenstate coupling matrices 

q(i) = R (fM q(i) R (q) ; ( 4 ) 

related to the corresponding weak-eigenstate coupling matrices, G\ q \ by means of the 
rotation matrices . The latters diagonalize the stop and sbottom mass matrices defined 
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in eqs.(8)-(9) of Ref.||. For the tt final states, we have 
6? [cot (3} = 



— im t (fi + At cot (3) 

im t (/i + A t cot (3) 



G 2 ^ [Caj S a , Sp\ 



-2Mi(T 3 W -gW^) s a+ ,+ 
(fis a + A t c a ) 



nit 

S/3 



+ Ac, 



^3 ^ [ C Q5 S «5 S /3] 



'2^ (if -Q^)c 



+ 



2m|s c 
S/3 



2* (-/iC Q + A t S Q ) 



(-/iC a + A t S a ) 

2rrt?s 



2M 2 z Q^s 2 w c a+/3 + 



(5) 



with s Q+/ 3 = sin(a + /3) etc. and Q( q \T$ the electric charge and 3rd component of 
weak isospin. For the b b final states, the following replacements are to be performed with 
respect to the G^[...] in eq.(EI): 



G? 
Gf 



G? [tan/?] , 
G-2 [s a , c a , —Cp] 
G3 ^ [s a , C a , Cp] . 



(6) 



The one-loop renormalized vertices for any of the relevant hadronic decays $ l — > qq are 
derived from the on-shell renormalized Lagrangian and can be parametrized in terms 
of two bare form factors K l L (q), K R (q) and the corresponding mass and wave-function 
renormalization counterterms Sm q and 8Z\ R associated to the external quark lines: 

gm q 



0\q) 



2M W 

the renormalized form factors being 



<(q) (l + 01(g)) P L + a R {q) (l + 0* R {q)) P R 



(7) 



O l L {q) 



(8) 



For each $* = A , h°,H° decaying into qq a straightforward calculation of the diagrams 
in Fig. lb yields a generic contribution of the form (summation is understood over a, b) 



Kl(q) 
K l R (q) 



871a s iCp 



r i ab 

(5) 



hra a iC F -4^r 



p(?) p(<?)* //^ n \ j_ p(?) T){<i)*ri _l m 9 0(9) 0(9)*/^ 

-"-2a -"-26 l°H ~~ U 12j + -n-la ^16 °12 H -"-la K 26 °0 



(9) 



The explicit expressions for the mass and wave-function renormalization counterterms are 
borrowed from Ref. || and will not be repeated here, and the various 3-point functions 
in eq.(H) have the arguments C... = C_(p,p', nig, m^, Triq b ) [fLOR ; they carry indices a, b 
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summed over. Finally, Cf = {N^ — l)/2Nc = 4/3 follows from summation over color 
indices. 

At the end of the day, the relative SUSY-QCD correction to each decay width of 
$' — > qq with respect to the corresponding tree-level width reads as follows, 

W = r{ %q)° iq) = + > ( 10 ) 

where T l (q) = — > qq) is the corrected width, and 

/ N C G F \ , . 2 _ . <,./■./„,.••>,. 
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(j = for i = 1 and j = 1 for % = 2, 3) , (11) 

are the corresponding tree-level widths. The upshot of our exhaustive numerical analysis 
is synthesized in Figs. 2-5. We treat the sbottom and stop mass matrices as in Ref. 

In Fig. 2a, where we fix M^o = 60GeV, we study the dependence of the SUSY-QCD 
correction ( |T0"D on the Higgs mixing mass parameter jj, for the three decays — > b b. We 
immediately gather that the sign of the correction is opposite to that of \i. For A and 
h° the correction is basically the same and can reach large values, e.g. \8§\ ~ 50% at 
tan/? = 30 and \p\ ~ 100 GeV. As stressed in Ref.g, the origin of the large SUSY-QCD 
contributions obtained in the presence of final states involving the 6-quark can be ascribed 



to their self-energy renormalization effects [11], which in our case go to the counterterm 



Snib/mb on eq.(f|). We remark that the corrections affecting the bb final states are larger 
for H° than for A and h°. The drawback, however, is that the huge effects obtained for 
H° — > bb at the highest values of tan (3 correspond to the smallest tree-level decay widths. 
In contrast, the less ambitious but still quite respectable quantum effects on A —> bb are 
larger the larger is its decay width. 

In Fig. 2b we study the alternative decays of A and H° into ti final states for parameter 
values in Region II. Here, the minimum value of the lightest stop mass has to be specified 
and we take mj = 65 GeV. Even though A — > bb is also dominant in Region II for the 
largest allowed values of tan {3 in this region, it has large QCD backgrounds. The heavy 
ti final states, however, are projected in the direction of the beam and can be identified 
through high-py leptons from semileptonic t-quark decays. Thus the heavy Higgs decays 
into t i final states, though they have a branching fraction smaller than that of the b b final 
states for tan/3 > 6, may in compensation be more manageable from the experimental 
point of view. For these decays we generally select more moderate values for tan f3 in 
order to make them sufficiently operative. From Fig. 2b we realize that of the two decays 
into ti, the most sensitive to SUSY-QCD radiative corrections is that of the CP-odd 
Higgs boson. Here, in contradistinction to the 6-quark final states, the main source of the 
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corrections lies in the structure of the form factors Kl,r on eqs.(||)-(^D - the top quark 
self-energies being negligible. 

Of course, we expect -and we have numerically verified- that the SUSY-QCD con- 
tributions drop off upon freely increasing the squark masses. However, in practice the 
asymptotic regime begins for fairly large values of these masses. For example, in Fig. 3 we 
study the Higgs decays into bb as a function of , for various tan/3. We see that, for 
tan/3 > 10, the corrections can reach several 10% even for m~ bi in the few hundred GeV 
range. 

Worth noticing is the asymptotic behaviour of the correction (|l(]) versus the gluino 
mass for the various Higgs decays. As shown in Figs.4a-4d, it takes a long time, so to 
speak, for the gluino to decouple. Corrections of ~ 50 — 60% for A°,h°,H° —>■ bb are 
obtained at high tan /3 from a mass value m- g ~ 150 GeV all the way out to 1 TeV - hence 
far beyond the present phenomenological bounds. In Figs.4c-4d we can also assess the 
dependence of A°,H° — > tt on m^, for fixed = 150 GeV and a moderate value of 
tan/3; and we see that even for stop masses as heavy as 100 GeV the corrections are longly 
sustained (above 10%) for practically any value of m g beyond the threshold singularities 
associated to points satisfying m g + m~ tl ~ m t . For gluino masses below these points, the 
corrections to A — > t i can be much larger. 

From Figs. 5a and 5b we read off the dependence of the SUSY-QCD corrections on M^o 
for different values of tan (3, and they are compared with the ordinary QCD corrections. 
We remind the reader that the QCD corrections to — > qq can be very large for light 
quarks |J . As in the decay of the charged Higgs [Q] , this is due to the appearance of 
a logarithmic term carrying a quark mass singularity, ~ log (M$</m g ), which stems from 
the anomalous dimension of the qq and qj^q operators. The complete one- loop (and 
renormalization group improved) formulae read as follows]] {b — 33 2n/ V 



r l (q) = r (q) 



1 /w \ 1 M * 

ba s (Mq>i) log 



2m q 



6?r 

4/bvr ( 



1 H (A$i + 3 log 



7r 2 m q 



(12) 

where the complicated functions are given by eqs.(3.7) and (2.26) of Ref.j| for 
i = l and i = 2,3 respect ivelyQ. From eqs. (|I2"D and (^T|) the standard QCD corrections 
3g = C^qcd(i) ~ Tq)/Fq to the various MSSM neutral Higgs decays can be computed and 
are included in Figs. 5a and 5b, where they can be compared with the SUSY-QCD effects 
(SI). 

From Fig. 5a we see that the decays — > b b receive negative standard QCD corrections 

1 This equation is equivalent to eq.(4.5) of Ref. except that we have corrected a missing factor of 2 
in the last logarithm. 

2 We have also corrected a missing factor of 2 in the third term on the RHS of eq.(2.27) of Ref.Q. 
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around 30 - 45% (Notice that we have plotted -tfj in Fig.5a.). For M A o > 100 GeV, 
the standard QCD correction to h° — > bb remains saturated at about —30% since the 
mass M^o also saturates at its maximum value, whereas the modes A°,H° —>■ bb obtain 
slowly increasing negative corrections. In contrast, A —>■ tt and H° — > tt, receive 
positive standard QCD corrections rapidly varying with the Higgs mass (Cf. Fig. 5b) 
Comparison with our Figs. 2-5 clearly shows that in many cases the supersymmetric effects 
are important since they can be of the same order as the standard QCD corrections. 

Overall, large SUSY-QCD quantum corrections are expected in the hadronic widths of 
the neutral Higgs bosons of the MSSM. They should be measurable, though with different 
techniques, both in e + e~ and in hadron machines. These supersymmetric effects can not 
only be comparable to the ordinary QCD corrections, but even dominant in some cases. 
Since they can have either sign, the net QCD correction would be found either much 
"larger" than expected, perhaps "missing" or even with the "wrong" sign; in any case, it 
should be revealing to hint at the SUSY nature of these higgses. Furthermore, we have 
found that, contrary to the situation with SUSY corrections on gauge boson observables, 
these effects decouple very slowly, especially with the gluino mass. Therefore, if SUSY is 
there, these corrections should also be there, and cannot be missed for a wide range of 
sparticle masses. However, Region II is out of reach of LEP 200, and even though part 
of the Higgs spectrum characterizing Region I is within its discovery range a complete 
experimental account of the MSSM Higgs sector will not be possible at LEP 200. In this 
respect, we have put forward the convenience of trying to see the kind of effects studied 
here in the large hadron machines, perhaps before an e + e~ supercollider be at work. In 
fact, the potentially large size of these effects indicates that they ought to be included 
in any serious analysis of supersymmetric Higgs production processes in hadron colliders. 
The combined information on branching ratios (from e + e~ ) and on cross-sections (from 
the Tevatron and/or at the LHC) should be very useful to pin down the nature of the 
observed effects. A more complete study should include the electroweak SUSY effects (in 
particular, the effect from the one- loop Higgs mass relations) , but as already mentioned in 
the companion paper || these are expected not to drastically alter the SUSY-QCD picture 
presented here. Our general conclusion is that quantum corrections on Higgs physics may 
be the clue to "virtual" Supersymmetry. 
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Figure Captions 

• Fig.l (a) Typical mechanisms for neutral Higgs production at hadron colliders; (b) 
SUSY-QCD Feynman diagrams, up to one-loop level, correcting the partial widths 
of A , h° } H° -> qq. 

• Fig. 2 (a) Dependence of the relative SUSY-QCD corrections S~(q) - eq.(|T0D - for 
i = 1,2,3 and q = b, upon the supersymmetric Higgs mass mixing term, /z, for light 
M A o = 60 GeV and given values of the other parameters (the scale of the abscissa 
is common to (b) below); (b) As before but for Sg' 3 (t) and heavy M40 = 400 GeV, 
for fixed m~ tl =65 GeV . 
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• Fig. 3 5~' 2 ' 3 (6) as a function of the lightest sbottom mass m~ bi for various tan/? and 
fixed ji and M A o. The other fixed parameters are as in Fig. 2. 

• Fig. 4 (a) Evolution of 5~(b) ~ ^|(^) (almost indistinguishable) in terms of the 
gluino mass for various tan/3. M40 is chosen light and the remaining inputs as in 
Fig. 2 (the scale of the abscissa is common to (b) below); (b) As before but for <5~(ft); 
(c) 5~(t) versus the gluino mass at fixed tan/3 and for various m~ ti (the scale of the 
abscissa is common to (d) below); (d) As before but for <5~(t). 

• Fig. 5 (a) 5~' 2,3 (&) for tan/3 = 30 (upper-born curves) and tan/3 = 4 (lower-born 
curves) as a function of M^o and the rest of the parameters as in Fig.2. The middle- 
born curves stand for the corresponding standard QCD corrections, 5 g , to the three 
decay processes. The latters being negative, we plot —5 g to ease comparison with 
the SUSY-QCD corrections at fi < 0; (b) As before but for <5~' 3 (t) and the range of 
M40 selected deep into Region II and three values of tan (3. The scale of the ordinate 
is common to (a). 
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